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Abstract We revisit a phenomenological description of turbulent thermal convec- 
tion along the lines proposed by Gough in which eddies grow solely by extract- 
ing energy from the unstably stratified mean state and are subsequently destroyed 
by internal shear instability. This work is part of an ongoing investigation for finding 
a procedure to calculate the turbulent fluxes of heat and momentum in the presence 
of a shearing background flow in stars. 

We describe the first steps towards a generalization of the time-dependent formu- 
lation of the mixing-length approach for radial pulsation, proposed by Gough [2], 
to nonradially pulsating stars. To this end, we incorporate a treatment of the influ- 
ence of a shearing background flow. In order to test and calibrate the formalism, 
we first compare its predictions with results of direct numerical simulation (DNS) 
of Rayleigh-Benard convection in air (Prandtl number, a = 0.71) with a strongly 
shearing background flow JTJ . We consider a plane-parallel layer of fluid confined 
between rigid horizontal perfectly conducting boundaries of infinite extent at fixed 
temperatures, the lower being hotter than the upper by AT. The upper boundary 
moves horizontally with constant velocity AU relative to the lower boundary, and at 
first we assume, in accordance with the Boussinesq approximation, that the shear, E, 
in the mean flow does not vary over the scale of an eddy. Previously we considered 
the local model only Q, which we extend here to account for nonlocal effects. 

Following J4] we solve the linearized dimensionless equations describing the 
dynamics in a statistically stationary flow of a viscous Boussinesq fluid confined 
between two horizontal planes. Linearized modes of convection are obtained by 
expansion about E := d z U = (d- = d/dz) to second order in E. The resulting ex- 
pressions for the eigenfunctions of the fluctuating temperature, T', and turbulent ve- 
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Fig. 1 Mean temperature profile T j AT compared with DNS data (left panel) and dimensionless 
convective heat flux (F c ) (right panel) as a function of height z for different AU values. 



locity field, m, = (u,v, w), are used to compute the turbulent fluxes of heat, F c °< wT', 
and momentum (Reynolds stresses), uju], in the manner of 0. A horizontal bar 
denotes the average over the horizontal plane. The local expressions are next aver- 
aged in the manner proposed by J3| to account for nonlocal effects: the eddies of 
vertical extent I centred at zo sample the temperature gradient over the spatial re- 
gion {zq — 1/2, zq + t/2), and the nonlocal fluxes, (F c ) and (HjUj), are determined 
from the contributions of all eddies centred between zo — 1/2 and zo + 1/2. These 
fluxes enter the conservation equations for heat, N = — d z T + (F c ), and momentum, 
S = d z U — (uw) / (J, which enable one to compute the mean temperature and velocity 
profiles, T(z) and U(z), across fluid layer. The total dimensionless heat flux (Nus- 
selt number), N, and total stress, S, are eigenvalues of the system. The diffusion 
approximation for the radiative flux is used. 

In Fig. Q] (left panel) we compare the mean temperature profile with the DNS 
data. Agreement of the functional form of the mean temperature profile between 
our model results and DNS data is obtained only for lesser shear in the model than 
in the numerical simulation. The right panel of Fig.[T]illustrates that, in agreement 
with the DNS data, the convective heat flux is reduced with increasing shear. We 
stress the need for further numerical simulations to enable us to make more detailed 
comparisons over a much broader range of Rayleigh and Prandtl numbers. 
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